Improvements  in  Code  TUNCOR  for 
Calculating  Wall  Interference  Corrections 
in  the  Transonic  Regime 

M.  H.  Rizk 
Flow  Industries,  Inc. 

21414-68th  Avenue  South 
Kent,  Washington  98032 


March  1986 


Final  Report  for  Period  April  1983  -  December  1985 


Approved  for  public  rainn:  distribution  unllmltad. 


ARNOLD  ENGINEERING  DEVELOPMENT  CENTER 
ARNOLD  AIR  FORCE  STATION,  TENNESSEE 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 


I 


NOTICES 


When  U.  S.  Govemnwnt  drawings,  specifications,  or  other  data  are  used  for  any  purpose  other  than  a 
definitely  related  Government  procurement  operation,  the  Government  thereby  incurs  no  responsibility 
nor  any  obligation  whatsoever,  and  the  fact  that  the  government  may  have  formuiated,  furnished,  or  in 
any  way  supplied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  regarded  by  impiication  or 
otherwise,  or  in  any  manner  licensing  the  holder  or  any  other  person  or  corporation,  or  conveying  any 
rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in  any  way  be  related 
thereto. 


Qualifled  users  may  obtain  copies  of  this  report  from  the  Defense  .Technical  Information  Center. 


References  to  named  commercial  products  in  this  report  are  not  to  be  considered  in  any  sense  as  an 
endorsement  of  the  product  by  the  United  States  Air  Force  or  the  Government. 


This  report  has  been  reviewed  by  the  Office  of  Public  Affairs  (PA)  and  is  releasable  to  the  National 
Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be  available  to  the  general  public,  including 
foreign  nations. 

APPROVAL  STATEMENT 


This  report  has  been  reviewed  and  approved. 


KEITH  L.  KUSHMAN 
Directorate  of  ITechnology 
Deputy  for  Operations 


Approved  for  publication; 

FOR  THE  COMMANDER 

LOWELL  C.  KEEL,  Lt  Colonel,  USAF 
Director  of  Technology 
Deputy  for  Operations 


MCU-HITv  CLASSIf  ICATIOW  DF  THIS  PAGE 


REPORT  DOCUMENTATION  PAGE 


llA  REPORT  SECURITV  CLASSIFICATION 

[lb  RESTRICTIVE  MARKINGS 

1  UNCLASS IFED 

_ 

z*.  SECURITY  CLASSIFICATION  AUTHORITY 

3.  PISTRIBUTION/AVAI lability  OF  REPORT 

Approved  for  public  release; 

distri- 

3b.  oeclassification/ooiwngrading  schedule 

button  unlimited. 

4  PEflFOnWING  ORGANIZATION  REPORT  NUMBEHIS) 

Is.  MONITORING  CHGANIZATION  REPORT  NUMeERfS) 

AEDC-TR-86-6 

6b.  name  of  performing  organization 

Flow  Industries ,  Inc . 

Sb.  OFFICE  SYMBOL 

tlfnppticobtct 

7b  name  op  monitoring  ORGANIZATION 

Be.  ADDRESS  rCily.  Statr  and  Xir  Ctidet 

21414  68th  Avenue  South 

Kent,  Washington  98032 

7b  ADDRESS  (Cfly.  SIbm  and  7IF  Codal 

NAME  OF  FUNOING/SPONSORING 

ORGANIZATION  Arnold  Engi¬ 
neering  Development  Centex 

B6.  OFFICE  SYMBOL 
(if  appLcabtet 

DOT 

9.  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

F40600-83-C-0003 

Be.  ADDRESS  fCify,  liNiN  and  ZIF  Codal 

10  SOURCE  OF  Funding  no$ 

Air  Force  Systems  Command 

Arnold  Air  Force  Station,  TN  37389 

PROGRAM 

E  LEMENT  NO 

PROJECT 

NO 

TASK 

NO 

WORK  UNIT 
NO 

1 1 1  TITLE  tlnelkdt  dacuntf  Clauiftealioni 

65807F 

Please  see  reverse  of  this  page.  I 

12  PERSONAL  AUTHORI5) 

Rizk.  M.H. .  Flow  Industries,  Inc.,  Kent.  Washington 


13a.  TYPE  OF  REPORT 

t3b.  TIME  COVERED 

14  DATS  OF  REPORT  fYt,  Mo..  Day/ 

16  PAGE  COUNT  1 

Final 

FROM  4/83  TO 

12/85 

March  1986 

35  1 

ie.  SUPPtEMENTA#«y  NOTATION 


Available  in  Defense  Technical  Information  Center  (DTIC). 


IT. 

COSATI  CODES  1 

18  SUBJECT  TERMS  <Cmi4fj)uc  on  rtrvane  tf  ntc€9»ary  and  idanitfy  by  bloci^  numter/ 

FIELD 

GROUP 

SUB  GR 

wind  tunnel  corrections  transonic  flow 

wall  interference 

20 

04 

1ft.  ABSTRACT  (CdnJMut  -or  rwMPi^  it  nmetMwry  and  Jdtfnft/y  by  Mi>cA  numberi 


Modifications  are  introduced  to  Code  TUNCOR  to  allow  its  use  in  de¬ 
termining  wall  interference  corrections  in  the  AEDC  IT  tunnel.  The  modif¬ 
ications  include  conversion  to  cylindrical  coordinates  and  converting  the 
measured  pressure  data  to  a  form  acceptable  by  the  code. 


20  QISTRIBUTION/AVAILABILITY  OF  ABSTRACT 

UNCLASSIFIED/UNLIMITED  □  SAME  AS  RPT  CS^DTlC  USERS  C 

21  abstract  SECURITY  CLASSIFICATION 

Unclassified 

22l.  NAME  OF  RESPONSIBLE  INDIVIDUAL 

22b  TZLePHONE  NUMBER 

22g  office  symbol 

(intiude  Area  Codei 

W.O.  Cole 

:615)  454-7813 

DOS 

EDITION  OF  1  JAN  IS  OBSOLETE 


IINr.lA?^RTi 


SECUniTV  CLASSIFICATION  OF  THIS  PAGE 


HCUniTV  CLJtSaiFKATION  OF  THIS  PAGE 


TITLE.  Concluded. 

Improvements  in  Code  TUNCOR  for  Calculating  Wall  Interference 
Corrections  In  the  Transonic  Regime 


niiViHiaiMu 


KCURITV  CLAUIPICATION  OF  THIS  FAQE 


PREFACE 


The  work  reported  herein  was  conducted  by  Flow  Industries,  Inc., 
21414  -  68tb  Avenue  South,  Kent,  Washington,  under  Air  Force  Contract 
F4060D-83-C-0003  with  Arnold  Engineering  Development  Center/PK,  Air 
Force  Systems  Command,  Arnold  Air  Force  Station,  Tennessee.  The  AEDC 
Project  Manager  for  this  work  was  Dr.  Keith  L.  Kushman,  AEDC/DOT. 

The  reproducibles  used  in  the  reproduction  of  this  report  were 
supplied  by  the  author. 


1 


AEDC-TR-Se-e 


CONTENTS 

Page 

l.O  INTRODUCTION  5 

Z.O  SUHHART  OF  MODIFICATIONS  TO  CODE  TUNCOR  7 

3.0  GOVERNING  EQUATIONS  10 

3. 1  Model  Governing  Equation  10 

3.2  Model  Boundary  Conditions  10 

4.0  FINITE  DIFFERENCE  EQUATIONS  22 

REFERENCES  26 

ILLUSTRATIONS 

Figure 

1.  Correction  Procedure  6 

2.  Geometrical  Configuration  Used  in  the  Original 

TUNCOR  Code  8 

3.  Cylindrical  Surface  on  Which  Pressure 

Measurements  are  Hade  3 

4.  Farallelopiped  on  Which  Body  Boundary 

Conditions  are  Applied  14 

5.  Body  and  Farallelopiped  Cross-Sectional 

Contours  15 

6.  Cylinder  on  Which  Body  Boundary  Conditions 

are  Applied  17 

7.  Body  and  Cylinder  Cross-Sectional  Contours  18 

8.  Cross  Section  of  Computational  Body  20 

APPENDIX 

A.  COMMENTS  ABOUT  COMPUTER  CODE  TUNCORC  27 

NOMENCLATURE 


3 


AEDC-TR-B6'6 


I.O  INTRODUCTION 

Wall  interference  effects  are  a  limiting  factor  on  the  accuracy  of  data 
obtained  from  transonic  vind  tunnel  tests.  Classical  linear  theory^  and 
modifications  to  this  theory  provide  an  insight  into  the  nature  of  these 
errors.  However,  sufficiently  accurate  estinates  for  practical  use  are  not 
obtainable  from  linear  theory.  Kemp^’^  formulated  a  procedure  that  uses 
nonlinear  transonic  computer  codes  to  determine  Mach  number  and  angle-of- 

g 

attack  corrections.  Murman  made  improvements  in  Kemp’s  formulation  and 
conducted  a  series  of  two-dimensional  airfoil  simulations  to  show  that  wall 
interference  corrections  are  posaible  for  strongly  supercritical  flows.  Rizk 
and  Murman  extended  the  approach  of  Kemp  to  practical  three-dimensional 
geometries.  They  developed  a  computer  code,  TUNCOR,  which  estimates  wind 
tunnel  wall  interference  corrections.  The  purpose  of  this  report  is  to 
describe  modificstions  to  code  TUNCOR  thst  increase  its  accuracy  and  allow  it 
to  be  compatible  with  the  AEDC  IT  wind  tunnel  measurement  system  geometry. 

In  addition  to  Mach  number  and  angle-of-attack  corrections,  the  transonic 
correction  procedure  used  in  code  TUNCOR  provides  an  estimate  of  the  accuracy 
of  the  corrections.  Lift,  pitching  moment  and  pressure  measurements  near  the 
tunnel  walls  are  required.  The  correction  procedure  may  be  divided  into  two 
steps.  In  the  first  step,  the  flow  about  the  test  model  is  simulated  numeri¬ 
cally  using  the  pressure  measurements  near  the  tunnel  walls  as  boundary 
conditions.  An  inviscid  transonic  flow  code  is  used  for  this  purpose.  In 
this  step  the  wing  and  tail  angles  of  attack,  cl.  and  ol.  ,  are  determined 

I  jW  Aft 

such  that  the  calculated  lift  and  pitching  moment  of  the  simulated  model  are 

equal  to  the  measured  lift  and  pitching  moment,  respectively.  The  angles  of 

attack  OL.  and  ol,  will  generally  be  different  from  the  experimental  values 
I  ,  w  T » t 

a  „  anda  .  This  is  due  to  the  viscous  effects  present  in  the  experiment  but  not 
e,w  fi ,  <■ 

in  the  numerical  simulation  of  the  flow  and  to  geometrical  differences  between 
the  test  model  and  the  simulated  model.  In  the  second  step,  the  flow  about 
the  model  in  free  air  is  numerically  simulated.  Angle-of-attack  corrections 
and  a  free-stream  Mach  number  correction  are  determined  such  that  the  calcu¬ 
lated  model  lift  and  pitching  moment  match  the  experimental  values  and  the 
calculated  Hach  number  difference  on  the  model  surface  in  the  tunnel  and  free 
air  is  minimized. 

A  summary  of  the  correction  procedure  is  given  in  Figure  1.  A  detailed 
description  of  the  procedure  is  given  in  Ref.  9> 
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Figure  1.  Correction  procedure 


6 


AEDC-TR-88-6 


2.0  SUMMARY  OF  MODIFICATIONS  TO  CODE  TUNCOR 

The  correction  procedure  requires  Che  numerical  simulation  of  the  flow 
about  the  test  model.  In  code  TUNCOR  this  simulation  is  achieved  through  the 
numerical  solution  of  Che  transonic  small  disturbance  equation.  The  formula¬ 
tion  upon  which  the  code  is  based  assumes  that  pressure  measurements  are  made 
on  the  boundaries  of  a  rectangular  parallelepiped  with  sides  close  to  the 
tunnel  walls  but  outside  the  boundary  layer  region  (see  Figure  2),  In  the 
first  step  of  the  correction  procedure,  Che  transonic  small  disturbance 
equation  is  therefore  solved  subject  to  the  boundary  conditions  applied  on  Che 
boundaries  of  the  rectangular  parallelepiped.  A  Cartesian  coordinate  system 
is  used. 

In  the  AEDC  IT  tunnel,  pressure  measurements  are  made  on  a  cylindrical 
surface  (see  Figure  3).  Conversion  to  cylindrical  coordinates  is  therefore 
necessary  to  allow  proper  numerical  simulation  of  the  tunnel  flow  subject  to 
the  measured  boundary  conditions. 

Code  TUNCOR  calculates  the  flow  in  the  tunnel  subject  to  Dirichlet 
boundary  conditions  for  the  perturbation  velocity  potential  ^  on  the  side 
boundaries.  It  is  therefore  necessary  to  convert  the  measured  data  into 
boundary  conditions  that  are  acceptable  by  the  code.  Interpolation  routines 
are  also  added  so  that  boundary  conditions  at  the  computational  mesh  boundary 
points  may  be  calculated. 

Finally,  improved  accuracy  is  obtained  by  replacing  the  transonic  small 
disturbance  governing  equation  and  the  first-order  boundary  conditions  by  the 
full  potential  equation  and  second-order  boundary  conditions,  respectively.  A 
detailed  description  of  the  governing  equation  and  boundary  conditions  is 
given  in  the  following  section. 
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3.0  GOVERNING  EQUATIONS 


3. L  FLOW  GOVERNING  EQUATION 

The  flow  is  governed  by  che  continuity  equation 

V'(pi)  -  0 

or  3  a  a 

gj-  (PUR)  +  tpv)  +  ^  (PwR)  =  0  (1) 

where  the  velocity  £  is  expressed  in  terms  of  the  perturbation  velocity  poten¬ 
tial  function  i|>  as  follows: 

a  =  V(X4<t») 

and  the  density  p  for  an  isentropic  flow  is  given  by 

p  =  [i  +  (1  -  (2) 

The  small  disturbance  assumption  made  in  Ref.  9  allows  expression  (2)  for  p 
to  be  expanded  in  terms  of  a  truncated  binomial  expansion.  When  this 
expansion  is  substituted  into  Eq.  (1),  the  transonic  small  disturbance 
equation  results.  The  full  expression  for  p  given  by  Eq.  (2)  is  used  in  the 
present  work. 

3.2  MODEL  BOUNDARY  CONDITIONS 

For  a  general  body  with  a  surface  defined  by 

B(X,Y,Z)  =•  0 

the  boundary  condition  specifying  zero  flow  through  the  body  surface  is  given 
by 

£*7B  =  0 

This  boundary  condition  may  be  expressed  as 

3y  3Y 

V  =  UR  ^  +  wR  ^  on  Y  -  Yg(X,2)  (3) 
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and 


V 


3Z  3Z 

_ s  V  _ s 

3X  R  3Y 


on  Z  =  Z^(X,Y) 


(4) 


vhere  the  surfaces 

V  =  V^(X,Z)  and  2  =  Z^(X,Y) 

specify  solid  boundaries* 

It  is  convenient  to  apply  boundary  conditions  on  coordinate  surfaces  when 
the  problem  is  solved  numerically.  For  surfaces  that  deviate  slightly  from 
the  coordinate  surface 

Y  =  Y 

o 

it  is  possible  to  make  the  following  Taylor  series  expansion: 


(P^)| 


=  (pv) 


lY  =  Y 


Y  =  Y 


-  -  »o>  ^ 


lY  =  Y. 


.  o[(Y.  - 


(5) 


Using  Eq.  (3)  and  Eq.  (1)  to  express  the  first  and  second  terms,  respectively, 
on  the  right-hand  side  of  Eq.  (5),  one  gets 


(Pv)J 


Y  =  Y 


=  3|[(WR)  (Y.  -Y„)] 

* 

+  0  j^{Y^  -  Y^)^j 


(pwR)  (Y  -  Y„) 
s  o 


] 


Y  =  Y 


For  surfaces  that  deviate  slightly  from  the  coordinate  surface 

2-2 

o 

it  is  possible  to  make  the  following  Taylor  series  expansion: 


(pwR) ; 


=  (pwR) 


IZ  =  Z 


z  -  z 


s 


+  0  1(2  -  Z 

I  s  o 


(6) 


(7) 
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Using  Eq.  (A)  and  Eq.  (1)  to  express  Che  first  and  second  terms,  respectively, 
on  the  right-hand  side  of  Eq.  (7),  one  gets 


(pwR)| 


(pUR)  (Z  -  Z  ) 
8  0 


]l^- 


W  [<»)  <2,  -  2„)]  Iz  . 

0  [(Z,-  z^)2] 


(8) 


The  usual  assumption  Chat  variations  in  the  Z  or  T  directions  are  negligible 
relative  to  variations  in  the  X  direction  is  made  here.  It  is  therefore 
assumed  that  Y  and  Z  derivatives  on  the  right-hand  sides  of  Eqs.  (6)  and  (8) 
are  Ol(Y^  -  Y^)^]  and  0[(Z^  -  Z^)^] ,  respectively.  The  boundary 
conditions  are  therefore 


(Pv)|^  = 


_i 

3X 


(PUR)  (Y 


(pwR) 


|z  =  Z^  ~  3X 


'  V 

s  o 


.  -  ’'«>] 
] 


lY  =  Y. 


(9) 


(PUR)  (Z,  -  Z) 
S  0 


-  Iz  - 


+  0 


['2.  - 


(10) 


The  usual  small  disturbance  boundary  conditions  are  derived  by  making  the 
following  Taylor  series  expansions: 


Y  =  Y  °  '^Iy  =  Y  *  ■  '’^o^ 

o  a 


(wR) 


iz  -  z  “ 


z  .  z  *  -  V 

8 


and  then  substituting  expressions  (3)  and  (4)  for  v  and  w  into  the  above 
expansions  after  setting  U  ■  1  and  neglecting  the  last  term  on  the  right-hand 
side  of  Eqs.  (3)  and  (4).  Therefore, 


and 


(11) 


(12) 


12 


AEDC-TR-Se-e 


Boundary  conditionB  (11)  and  (12)  have  been  used  in  Ref.  9;  however,  the  more 
accurate  conditions  (9)  and  (10)  are  used  in  this  work. 

In  Ref.  9,  a  Cartesian  computational  mesh  is  used  in  the  calculations. 

The  small  disturbance  boundary  conditions  for  the  wing  and  tail  are  applied  on 
planar  mean  surfaces  in  the  usual  manner.  The  body  boundary  conditions  are 
also  applied  on  the  planar  surfaces  of  a  rectangular  parallelepiped  (see 
Figure  4).  The  boundary  conditions  at  these  surfaces  are  given  by 

3(t)P(x.e)  _  3<)\x.9)  3sNx.e)/aB  , 

^  asP(e)/36 

where  S^(x,6)  is  Che  path  length  along  the  body  cross-sectional  contour  at 
a  given  x  station,  5^(9)  is  the  corresponding  path  length  along  the  paral¬ 
lelepiped  cross-sectional  contour,  n  is  the  coordinate  in  the  direction  normal 
to  the  body  contour,  and  N  is  the  coordinate  normal  to  the  parallelepiped 
contour  (see  Figure  5).  Relation  (13)  between  the  boundary  condition  applied 
on  the  rectangular  parallelepiped  and  chat  applied  on  the  body  surface  is 
chosen  to  match  the  mass  flux  introduced  by  the  two  conditions.  Let  the  body 
cross-sectional  contour  at  station  x  be  defined  by 

r  -  r’'(x,e) 


The  value  of  3s'*(x,9)/30  can  then  be  expressed  as 


3S^(x, 9) 

b 

cos  V 

where 


and  T  is  the  angle  between  the  body  contour  and  the  radial  direction  as 
indicated  in  Figure  5.  Similarly,  if  the  parallelepiped  cross-sectional 
contour  is  defined  by 

r  =  rP(e) 


then  the  value  of  3S^(6)/36  is  expressed  as 


3SP(0) 

36 


cos 
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Cross  Section 


Figure  5.  Body  and  paralleloplped  cross-sectional  contours 
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where 

=f  -  tP 


and  is  Che  angle  between  Che  parallelepiped  contour  and  the  radial 
direction  as  indicated  in  Figure  5>  Boundary  condition  (13)  applied  at  the 
parallelepiped  surface  may  Chen  be  rewritten  as 


a4(x,0,rP) 

3N 


3<t(x,P.r^)  r^ 


cos 


3n 


COB 


<U) 


As  indicated  in  Figure  4,  Che  rectangular  parallelepiped  extends  over  the 
length  of  the  computational  domain.  Boundary  condition  (14)  is  valid,  however, 
only  in  Che  interval  along  which  the  body  extends  (i.e.,  between  Sections  B  and 
C).  Beyond  this  interval,  Che  right-hand  side  of  Eq.  (14)  is  replaced  by  zero. 

The  corresponding  boundary  condition  formulation  in  cylindrical 
coordinates  replaces  Che  rectangular  parallelepiped  by  a  cylinder  (see 
Figure  6).  Eq.  (14)  is  applied  in  Che  body  region  between  Sections  B  and  C. 
Here  the  boundary  condition  reduces  to 


- K - *-5^“ — 


(15) 


r  cos  V 
o 


where  is  the  radius  of  the  cylinder  (see  Figure  7).  If  the  body  is 
nearly  axisyrsmetric,  its  cross-sectional  contours  will  deviate  only  slightly 
from  circular  cross  sections  and,  therefore,  cos  will  be  approximately 
equal  to  unity.  In  this  case,  cos  may  be  replaced  by  1  in  Eq.  (15). 
Upstream  of  Section  B  and  downstream  of  Section  C  the  condition 


3(|i(x,6,r  ) 

is  applied. 

The  application  of  the  body  boundary  conditions,  in  the  manner  described 
above,  on  a  parallelepiped  or  a  cylinder  modifies  the  problem  that  was 
originally  to  he  calculated.  This  formulation,  however,  does  simplify  the 
numerical  calculations  and  is  to  be  viewed  as  a  crude  approximation  allowing 
the  introduction  of  body  effects  into  the  problem  without  much  complexity. 
Refinements  that  better  approximate  the  original  problem  are  possible  and  will 
be  briefly  described  below.  First,  however,  a  description  of  the  modified 
model  that  corresponds  to  the  above  formulation  is  given. 
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The  application  of  the  boundary  condition 

9i|>(x,6,r  ) 

upstream  of  Section  B  and  downstream  of  Section  C  corresponds  to  the  intro¬ 
duction  of  a  cylinder  in  these  regions.  Therefore,  the  flow  cannot  penetrate 
the  surface  of  the  cylinder.  This  is  in  contrast  to  the  original  problem  in 
which  the  flow  is  not  restricted  upstream  or  downstream  of  the  body.  In  the 
body  region,  between  Sections  B  and  C,  the  new  formulation  corresponds  to  a 
change  in  Che  original  body  shape.  In  the  new  formulation,  the  body  becomes 
wrapped  around  the  cylinder  as  indicated  in  Figure  8.  This  would  indicate 
that  as  the  cylinder  radius  reduces  its  effects  in  the  body  region  and  outside 
the  body  region  also  decrease. 

Assuming  that  Che  body  defined  by 

r  “  r^(x,6) 


is  nearly  axisymmetric ,  then  Che  formulation  chat  includes  the  cylinder  of 
radius  r^  will  modify  the  boundary  in  Che  body  region  so  Chat  it  is  given  by 

r  =  r^(x,e)  = 

This  formula  indicates  Chst,  as  r^  approaches  0,  the  modified  body  approaches 
Che  original  body.  Similarly,  in  the  case  of  Cartesian  coordinates,  the  modi¬ 
fied  body  boundary  is  given  by 


Although  not  implemented  in  code  TUNCOR,  the  following  suggestions  which 

i 

give  a  better  approximation  to  Che  original  problem  are  briefly  discussed 
here.  In  the  case  of  Cartesian  coordinates,  the  rectangular  parallelepiped 
may  be  eliminated  upstream  of  Section  B  and  downstream  of  Section  C.  In  this 
case,  one  will  apply  the  body  boundary  conditions  on  the  sides  of  the 
parallelepiped  in  the  body  region.  This  will  better  simulate  the  original 
body.  Another  suggestion  is  to  eliminate  the  rectangular  parallelepiped 
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upstream  of  Section  B  only.  In  this  case,  boundary  conditions  can  be  imposed 
downstream  of  Section  C  to  model  Che  wake  region.  The  same  suggestions  apply 
to  cylindrical  coordinates;  however,  an  additional  complication  occurs  in  this 
case.  The  governing  equation  in  cylindrical  coordinates  is  singular  at  the 
axis  and,  therefore,  cannot  be  used  to  find  Che  solution  there.  In  this  case, 
one  may  replace  Che  cylinder  outside  Che  body  region  by  a  cylinder  with  a 
small  diameter  to  reduce  its  effects  there,  or  else  one  would  be  required  to 
solve  Che  governing  equation  on  Che  axis  in  Cartesian  coordinates  while 
solving  it  elsewhere  in  cylindrical  coordinates.  The  modifications  suggested 
above  will  increase  the  accuracy  of  simulating  Che  body;  however,  they  would 
also  increase  Che  complexity  of  solving  Che  equations. 
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4.0  FINITE  DIFFERENCE  EQUATIONS 


The  aolution  of  the  governing  equation  subject  to  the  specified  boundary 
conditions  is  found  numerically.  Therefore,  a  mesh  of  discrete  points 
with  mesh  spacings  AX,  AY,  AZ  is  defined  in  the 
computational  domain  with 

1  £  i  ^  I 

1  1  j  1  J 

1  <  k  <  K 


All  Neumann  boundary  conditions  are  applied  on  half-mesh  coordinate  surfaces 
(i.e..  Y  “  or  Z  =  We  shall  refer  to  mesh  points 

neighboring  boundary  surfaces  by  boundary  neighbor  mesh  points.  Other  mesh 
points  away  from  boundaries  are  referred  to  as  interior  mesh  points. 

Discretizing  Eq.  (1)  at  the  interior  mesh  point  (i,j,k)  leads  to  the 
following  finite  difference  equation: 


+ 


+ 


-  (pUR)' 


AX 

AY 

AZ 


-  0 


This  equation  may  be  rewritten  as 

a’’ 

5x 


a’’  ■  •  i.Jl 

+  (pwRj^'J’*^  2 


(16) 


3^  3^  3^ 

where  denote  backward  difference  quotients. 

The  finite  difference  approximation  to  the  governing  equation  at  a  mesh 
point  (ijjik)  neighboring  the  computational  boundary 


Y  -  Y 

o 
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is  given  by  Eq.  (16).  However,  now  in  this  equation  is 

replaced  by  the  finite  difference  approximation  to  Eq.  (9),  where 

Y  =  ^  AY 

o  —  z 


Similarly,  the  finite  difference  approximation  to  the  governing  equation  at  a 
mesh  point  (i,j,k)  neighboring  the  computational  boundary 

Z  »  Z 

o 


is  given  by  Eq  .  (16).  However,  now  ( (X#R)^ ^  in  this  equation  is 
replaced  by  the  finite  difference  approximation  to  Eq.  (10),  where 


Z  = 
o 


The  finite  difference  approximation  to  the  governing  equation  at  a  general 
mesh  point  (i,3,fc)  may  be  written  as 


.j.k 


*|^(PV 


■  0 


(17) 


where  the  value  of  » j is  zero  unless  the  mesh  point  (X^  ^  ,Y^ , Z*^) 

neighbors  the  computational  boundary 


Y  =  Y 


in  which  case 


1  ”1 

8  O 

s  V  ^  A 

AY 

i  i  V  1 

AY 

■  ^  ’^o  ^  2 

=  j 

Y  -  Y 
s  o 

Y^  =  Y  -  i 

AY 

( 

AY 

'  ^  ^0  2 

neighbors  the  computational 

boundary 

Z  =  2 

o 
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in  which  case 


Z  -  Z 
S  Q 

AZ 


=  2  * 
o 


1 

2 


tz 


The  values  of  and  are  1  except  on  computational  boundaries  normal  to  the 
y  and  Z  coordinates,  respectively,  where  they  are  zero. 

Finally,  writing  Eq.  (17)  in  terms  of  the  perturbation  potential  with 

U  =  1  + 

^  -  R  ♦y 
„  = 

and  making  the  transformations 

K  “  W 

C  "  CCY) 
r  -  ri(2) 


with  the  objective  of  using  a  stretched  mesh  away  from  the  aerodynamic  model 
(at  the  model  C  =  X,  C  -  V,  ri  =  Z) ,  the  following  is  the  resulting 
equation 


T  .  1  . 

)1 


<1  + 


=  0 


(18) 


.  a*^ 

where  ^  denote  centered  difference  quotients  and  p  is  the  density 

modified  by  the  artificial  viscosity  required  for  the  stability  of  the 
difference  scheme. 
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The  above  modifications  have  been  introduced  into  computer  code  TUNCOR 
(see  the  appendix),  which  predicts  Kach  number  and  angle-of-attack  corrections 
for  wind  tunnel  models.  The  code  has  been  submitted  to  AEDC  where  it  is 
presently  being  tested. 
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Code  TUNCOR  <Ref.  9)  has  been  developed  to  evaluate  wind  tunnel  wall 
interference  corrections  for  three-dimensional  aircraft  (wing-body-tail) 
configurations.  The  modified  version  of  this  code  (TUNCORC)  allows  the  option 
of  using  a  cylindrical  coordinate  system  in  addition  to  the  original  Cartesian 
coordinate  option.  Also,  code  TUNCORC  allows  the  use  of  the  full  potential 
equation  with  second-order  boundary  conditions  in  addition  to  the  original 
option  of  using  the  small  disturbance  equation  with  first-order  boundary 
conditions.  Finally ,  code  TUNCORC  includes  subroutines  that  transform  the 
measured  pressure  data  (on  either  a  cylindrical  surface  or  planar  surfaces)  to 
the  corresponding  potential  function  at  the  computational  boundary  for  use  in 
solving  the  tunnel  boundary  value  problem.  The  code  evaluates  Mach  number  and 
angle-of-attack  corrections  by  the  theory  developed  in  Ref.  9  and  in  the  main 
text  of  this  report.  Code  TUNCORC  includes  a  user's  manual.  The  comments 
given  here  are  not  intended  to  give  a  complete  description  of  the  code; 
however,  they  complement  the  user's  manual  included  in  code  TUNCORC. 

Code  TUNCORC  uses  a  simple  Cartesian  (or  cylindrical)  computational  mesh. 
The  mesh  point  coordinates  are  (X^,  Y-^,  Z^)  where 

1^1,2,.  .  .  |X 
j-1,  2]  .  .  .  ,J 
k  =  1,  2,  .  .  .  ,  X 

The  computational  mesh  is  uniform  inside  a  rectangular  parallelepiped  (or 
cylinder)  enclosing  the  model.  Outside  the  parallelepiped  (or  cylinder),  the 
mesh  is  stretched.  The  model  surface  boundary  conditions  are  applied  on  half¬ 
mesh  planar  surfaces  (i.e.,  X  =  or  Y  =  or  2  ■  2*^*^^^^^). 

This  requires  a  number  of  modifications  to  the  input  geometrical  model  config¬ 
uration.  In  other  words,  due  to  the  simplicity  of  the  computational  mesh,  the 
flow  is  generally  computed  about  a  modified  model  that  differs  slightly  from 
the  input  model.  A  discussion  of  the  nature  of  these  modifications  is  given 
in  this  appendix. 

The  user  has  the  freedom  to  specify  the  mesh  spacings  AY  and  AZ.  The  mesh 
spacing  Ax,  however,  cannot  be  arbitrarily  specified  if  the  sweep  angle  of 
the  wing  leading  edge  is  nonzero.  The  mesh  spacing  AX  is  restricted  to 

.  AZ 

Ax  =  tan  lb  ,  n  -  1 ,  2 ,  .  .  . 
n  w 
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This  restriction  on  fix  allows  the  wing  leading  edge  to  pass  through  the  mid¬ 
point  between  neighboring  mesh  points  in  the  X  direction  as  shown  in  Figure  A-1. 
In  the  case  of  a  leading-edge  sweep  angle  of  zero,  the  user  may  specify  fix. 

Once  fix  is  determined,  the  tail  leading  edge  is  restricted  to  one  of  the 
following  sweep  angles. 


Therefore,  the  code  will  adjust  the  input  sweep  angle  6^  to  the  closest 
allowable  value  i|/^.  The  tail  is  also  displaced  in  the  Z  direction  so  that  its 
base  lies  on  the  parallelepiped  (or  cylindrical)  surface  on  which  the  body 
boundary  condition  is  applied  (see  Figure  A-1).  This  parallelepiped  is  speci¬ 
fied  by  the  user;  however,  the  parallelepiped  also  is  modified  so  that  its 
planes  lie  at  the  mid-distance  between  two  neighboring  mesh  planes. 

As  part  of  the  output,  the  code  prints  any  adjusted  value  and  the  corres¬ 
ponding  input  value. 

The  input  data  are  shown  below  for  the  calculated  example.  In  this 
example,  all  dimensions  are  evaluated  in  inches.  The  user,  however,  has  the 
freedom  of  choosing  the  length  units.  Pressure  measurements  in  the  wind 
tunnel  are  given  in  the  input  data  in  terms  of  the  corresponding  pressure 
^®®^ficientB.  The  measured  lift  XLFTE  and  the  measured  moment  XMOME  are 
normalized  by  the  total  pressure,  so  that  XLFTE  and  XMOME  have  the  dimensions 
(unit  length)  and  (unit  length)  ,  respectively.  The  model  lift 
normalized  by  the  total  pressure  is  given  by 


L 

P 

o 


/  - 

model  (1  + 
surface 


n»k  dS 


while  the  model  lift  normalized  by  is  given  by 


*  *2 
Q  U 

^OD  OP 


f 

YHi  J 


n*k  dS 


model 

surface 


(:  .  ^ 


The  first  expression  is  used  in  the  code  since  it  is  independent  of  free-stream 
conditions  (M^, 
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HOMENCLATURE 


S.d<n*>  - 


k 

L 


e,d 


da 


“K 


a.t 


a,w 


e,d 

/<>v  -  v’ 

/H^d, 

Ur.t'EiW 

■  unit  vector  vertically  upward 
=  lift/P_ 


=  lift 

=  pitching  moment  about  the  axis  x  =  x  /P 

m  o 

=  Mach  number 

=  unit  vector  normal  to  model  surface 

■  *'V'»’  “f,»‘  “f.t* 

^  “  total  pressure 

=  velocity  normalized  by  free-stream  velocity  magnitude 

_  U  ,  Cartesian  coordinates 

|r  ,  cylindrical  coordinates 

v,w  “  perturbation  velocity  components  normalized  by  free-stream 
velocity  magnitude  in  X,Y,Z  directions 

=  1  +  u 

**  velocity  component  in  the  X  direction 
x,y,z  =  Cartesian  coordinate  system 
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X 

Y 

Z 

a 

Y 

Aa 

P 

P 

<l> 

* 


cylindrical  coordinate  system 

X  coordinate  of  axis  about  which  the  model  pitching  moment  is 

X 

I  y  ,  Cartesian  coordinates 

I  6  ,  cylindrical  coordinates 

I  z  ,  Cartesian  coordinates 

‘r  ,  cylindrical  coordiates 

angle  of  attack 

ratio  of  specific  heats 

Mach  number  correction 

angle-of-attack  correction 

density  normalized  by  free-stream  density 

density 

perturbation  velocity  potential  for  the  tunnel  flow 
perturbation  velocity  potential  for  the  free-air  flow 


Subscripts 


f 

F 

H 

s 

t 

T 


model 

measured  quantity  or  experimental  condition 
corrected  condition 

calculated  quantity  for  the  model  in  an  inviscid  free-air  flo( 
Mach  number 
model  surface 
tail 

calculated  quantity  for  the  model  in  an  inviscid  tunnel  flow 
wing 

undisturbed  condition 


measured 


31 


